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In recent years there has been growing interest in verifying the horizon-scale homogeneity of the 
Universe that follows from applying the Copernican Principle to the observed isotropy. This program 
has been stimulated by the discovery that a very large void, centred near us, can explain supernova 
luminosity distance measurements without dark energy. It is crucial to confront such models with 
as wide a variety of data as possible. With this application in mind, I develop the relativistic theory 
of linear scalar perturbations on spherically symmetric dust (Lemaitre-Tolman-Bondi) spacetimes, 
using the covariant 1 -f 1 -f 2 formalism. I show that the evolution of perturbations is determined by 
a small set of new linear transfer functions. If decaying modes are ignored (to be consistent with the 
standard inflationary paradigm), the standard techniques of perturbation theory on homogeneous 
backgrounds, such as harmonic expansion, can be applied, and results closely paralleling those of 
familiar cosmological perturbation theory can be obtained. 

PACS numbers: 98.80.Cq, 98.80.Jk 
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I. INTRODUCTION 

Since shortly after the discovery of the cosmic mi- 
crowave background (CMB) radiation, it has been clear 
that the Universe is very nearly isotropic, with depar- 
tures from uniformity in the CMB temperature at the 
level of 1 part in 10^ (apart from the dipole of presum- 
ably kinetic origin). Since the CMB comes to us from the 
greatest visible distances, that observed isotropy tells us 
about the symmetry of spacetime on the largest observ- 
able scales. The natural assumption is that this isotropy 
cannot be an accident of our location, and that any ob- 
server would see essentially the same picture, regardless 
of where they were located. This notion is known as the 
Copernican Principle, and implies that the Universe is 
essentially homogeneous on the largest observable scales. 

However, cosmologists are in the business of determin- 
ing the nature of the Universe through observation and 
deduction, rather than through philosophical postulate. 
As unappealing as it may appear, a Universe with sub- 
stantial radial inhomogeneity on the largest scales is con- 
sistent with the observed isotropy. Galaxy surveys are a 
means to quantify homogeneity, but are limited in their 
reach (see, e.g., [1, 2]), and radial inhomogeneity is dif- 
ficult to disentangle from redshift-dependent effects such 
as evolution. Direct observational constraints on homo- 
geneity over distances comparable to that to the last scat- 
tering surface, which emitted the CMB, are unavailable. 

A few novel proposals exist in the literature for obser- 
vational signatures of radial inhomogeneity. The CMB 
radiation we observe is partly scattered from within our 
past light cone, and this may allow us to constrain in- 
homogeneity [3, 4]. A consistency relation between lu- 
minosity distances and the Hubble rate, which must 
be satisfied for homogeneous and isotropic Friedmann- 
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Robertson- Walker (FRW) models, has been discussed [5]. 
Recently the time drift of cosmological redshifts was pro- 
posed as a test of homogeneity [6] . 

This might appear as something of a "dotting the i's 
and crossing the t's" cosmological exercise, if not for the 
realization nearly ten years ago that a particular form 
of radial inhomogeneity can actually mimic the effect of 
accelerated expansion in supernova data [7] . The model 
studied was the spherically symmetric dust, or Lemaitre- 
Tolman-Bondi (LTB) spacetime [8, 9, 10]. Since that 
time many studies have confirmed and elaborated upon 
this idea (see, e.g., [11, 12, 13, 14]; more references can 
be found in the brief review [15]). The basic idea is quite 
simple: the increasing expansion rate in time that stan- 
dard dark energy models describe is replaced with in- 
creasing expansion rate towards the centre of a radial 
inhomogeneity. These two possibilities are difficult to 
distinguish if our observations are limited to the surface 
of our past light cone. 

To get the sign of the effect right (apparent accelera- 
tion rather than deceleration), we must be near the centre 
of an underdensity, or void. The void must be nonlinear 
today if it is to mimic the zeroth-order effect of accel- 
eration. Additionally, the void must extend to a great 
distance, namely several hundred Mpc, corresponding to 
the redshifts at which the relevant supernova data lies, 
z ~ 0.2 — 1. The presence of such a large void is strongly 
at odds with the standard picture of structure formation, 
where the largest nonlinear voids today appear at scales 
of tens of Mpc. Nevertheless, the difficulties with under- 
standing the theoretical basis for dark energy, as well as 
the coincidence problem, have led many researchers to 
examine such void models as alternatives to dark energy. 

In the context of void models for acceleration, and, 
more generally, considering the more fundamental is- 
sue of large-scale homogeneity, it is very important to 
constrain radial homogeneity using actual observations. 
Clearly, the more data we use to do this the better. LTB 
models contain one (technically two; see below) free ra- 



dial functions which could be adjusted to optimize the 
fit to any single data set. Therefore it is with comple- 
mentary data that our greatest hope lies in constraining 
such models. Surprizingly, few studies of void models 
for acceleration have gone beyond fitting to supernova 
data. Notable exceptions are Refs. [11, 13], who consid- 
ered CMB data, and, very recently, Rcf. [14], who con- 
sidered a wider range of data including CMB and baryon 
acoustic oscillations. 

One important class of data which has remained out of 
reach so far involves the behaviour of scalar density per- 
turbations, i.e. the structure we observe on all scales. The 
evolution of this structure at linear order is necessarily 
affected by the presence of a large radial inhomogeneity 
or void in the zeroth-order background. While luminosity 
and angular diameter distances and redshifts are simple 
to calculate in LTB spacetimes, no one has yet studied 
perturbed LTB spacetimes. It may indeed appear daunt- 
ing to tackle perturbations on a background spacctimc 
which is not itself homogeneous [40]. However, with the 
appropriate tools, the task turns out to be not too much 
more difficult than the familiar perturbation theory on 
FRW backgrounds. 

Here I develop a formalism to study the evolution of 
linear scalar perturbations on LTB backgrounds. The 
approach is fully general relativistic, employing the co- 
variant 1-1-1-1-2 formalism (sec, e.g., [16, 17, 18]), which 
generalizes the standard 14-3 approach (see [19] for a 
recent review). This approach, while less familiar than 
the metric-based approach to many, offers several advan- 
tages, both in providing concise dynamical expressions 
which manifestly respect the symmetries of the problem, 
as well as in emphasizing physical and in principle ob- 
servable quantities. 

I do not develop the theory of perturbations on LTB 
backgrounds in complete generality. Rather, I target the 
development towards application to realistic void profiles 
that are able to mimic acceleration. My only approxima- 
tion is to ignore the coupling between scalar and tensor 
modes that curved backgrounds generally facilitate. I 
argue that for the applications of interest, the sourccd 
tensors will have a negligible effect on the scalar dynam- 
ics. I determine only the 2-scalar degrees of freedom, but 
these should be the most relevant observationally. While 
the background is taken to be spherically symmetric, the 
perturbations are completely free (apart from the linear 
and scalar approximations). 

It is important to choose void profiles that can fit 
within the standard inflationary paradigm, according to 
which at early times the Universe is expected to be very 
homogeneous. [One possible origin for a large void (and 
perhaps the most conservative) is as a rare, large am- 
plitude fluctuation within the primordial Gaussian ran- 
dom field produced during inflation. The only other ori- 
gin discussed for such a void appears to rely on highly 
speculative quantum cosmological ideas [20]]. Therefore, 
such a void cannot contain a significant decaying mode 
today, without spoiling primordial homogeneity. It is dif- 



ficult enough to believe that such a large void exists and 
is centred near us, without compounding the difficulties 
by giving up the standard inflationary scenario as well! 
Dropping the decaying mode also reduces our freedom 
to tweak an LTB model to fit data, so this should yield 
more predictive models. 

The dropping of the decaying mode provides a tremen- 
dous advantage in treating perturbations. Namely, since 
at early times the spacetime approaches FRW with linear 
fluctuations (one of which is the large proto-void itself) , 
all of the familiar tools used in perturbation theory on 
FRW can be applied at those early times, such as har- 
monic expansion and the use of statistical relations satis- 
fled by primordial Gaussian random flelds. At late times, 
when the LTB background has become nonlinear, it will 
be possible to relate the scalar fluctuations to the initial 
conditions using a small set of new LTB linear trans- 
fer functions. It should be possible to address many of 
the important observational signatures of structure using 
this formalism. The current work lays out the formalism; 
applications to speciflc observations will follow [21]. 

I begin in Sec. II with a brief review of the relevant co- 
variant 1-1-3 formalism. Next I describe the application 
of the further 1 + 1-1-2 decomposition to LTB models, 
before describing the exact LTB solution. Section III lays 
out the perturbation formalism, beginning with a careful 
definition of perturbations, before introducing the LTB 
transfer functions, and then setting up the initial condi- 
tions. Finally, Sec. IV presents some numerical results. I 
use signature (— , +, +, +) and set c = 1 throughout. 



II. 1 -M + 2 DECOMPOSITION FOR LTB 
BACKGROUND 



In the 1-1-3 covariant approach to relativistic cosmol- 
ogy (see, e.g., [19] for a recent review), all quantities are 
decomposed with respect to a fundamental congruence 
of timelike worldlines, with tangent vector field u''. A 
spacetime containing a single type of matter (pressure- 
less dust in the case of LTB) is perfectly suited to a co- 
variant treatment using the 1 + 3 decomposition, since 
the comoving worldlines naturally define a special time- 
like congruence, at least until worldlines intersect. When 
the congruence is twist-free, it is hypersurfacc-orthogonal 
[22] , and hence we can define spacelike slices everywhere 
orthogonal to the congruence. 

The spherical symmetry of the LTB background is well 
suited to further decomposing matter and metric quan- 
tities according to a 1 + 1 + 2 decomposition, since a 
spacelike congruence on each slice is naturally provided 
by the radial direction. In this section I provide a covari- 
ant 1 + 1 + 2 treatment of the LTB background, along 
the lines of that presented in Refs. [17, 18, 23]. This will 
set the stage for the treatment of linear perturbations in 
Sec. HI. 



A. 1 + 3 decomposition 

To begin, we will need a set of quantities which com- 
pletely describe the spacetime in terms of the 1 + 3 de- 
composition. Henceforth u^ will always be taken to be 
the timelike vector field tangent to the comoving world- 
lines, with normalization m^'m^j = — 1. The tensor 



/i^, = (5"„ + u^u. 



(1) 



projects orthogonal to w^. When the twist vanishes, we 
can consider /i^^ to be the metric tensor for the orthog- 
onal spacclike slices, and we can also use hf^, to define a 
spatially projected covariant derivative according to 
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for any tensor T^,-^^^...^^ orthogonal to u^ in all of its 
indices. To describe the time evolution of any covariant 
quantity, we employ the proper time derivative along the 
comoving worldlines. 



We will also need two tensors derived from the Weyl 
tensor, C^^Xp, which characterizes the nonlocal part of 
the gravitational field. The electric, £'^,y, and magnetic, 
H^^, parts of the Weyl tensor are defined by 
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The electric Weyl tensor describes the nonlocal tidal 
gravitational field, while the magnetic part describes, 
at least at linear order, propagating gravitational waves. 
Both Efj^i, and H^^i, are fully orthogonal to w^. For the 
case of a homogeneous and isotropic FRW cosmology, 
we have 9 = 3H, where H is the Hubble rate, and 
c^!/ = '-^pi' = fl/j = Efj^i, ~ Hf^y = 0. 

If the pressure (and hence the anisotropic stress) van- 
ishes, then the matter content is described by the simple 
stress-energy tensor 
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for any tensor quantity Tf^,^...p. Finally, angled brack- 
ets around tensor indices indicate the spatially projected, 
symmetric, and trace-free part: 



and the tensorial curl is defined by 



(4) 



with p the local energy density as judged by comoving ob- 
servers. The spatially projected stress-energy conserva- 
tion law, h'^jjT^^'^ ~ 0, i.e. the momentum conservation 
law, becomes 



a^ = 



(9) 



(as long as p ^ 0), representing the familiar result that 
comoving dust worldlines arc geodesic. 

The above kinematical, gravitational, and matter dc- 
(5) grces of freedom provide a complete description of the 
spacetime. It only remains to specify their evolution. 
Here e^^\ = e^^x^.u'^, where e^„\^ is the totally antisym- The twist ujf^^ of the congruence is a special case, and 



curir,,^ = e 
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metric four- volume element. 

We can characterize the comoving congruence kine- 
matically via the following standard decomposition of the 
covariant derivative u,, 



satisfies the evolution equation [19] 
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where 



e/ji/A'j-''''*'/2. Therefore, if the twist vanishes 



The scalar 6 measures the local volume rate of expansion 
of the congruence, while trace-free, symmetric tensor af^^^ 
and antisymmetric u;^^ measure the local rates of shear 
and twist of the congruence, respectively. The vector 
a^ = u^-^puP measures the acceleration of the comoving 
worldlines. Each of a^^, oj^,,^, and a^ are orthogonal to 
m'' in all of their indices. 



initially, it vanishes for all time. Henceforth the twist will 
always be set to zero. This amounts to ignoring the pres- 
ence of vector perturbation modes at linear order, which 
is a reasonable assumption within the standard inflation- 
ary scenario. Also, this means we can define spacelike 
slices orthogonal to the comoving flow. 

The evolution of the remaining quantities is deter- 
mined by the following set of equations [19]: 



Energy conservation: 

Ray chaudhuri : 

Shear evolution: 



Efj^i, evolution: 



Hfj,i, evolution: 



P = 
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(11) 
(12) 
(13) 
(14) 
(15) 



The evolution of each quantity is determined by a damp- 
ing term, proportional to the expansion (possibly cor- 
rected by the shear), together with coupling terms, such 
as the direct coupling between shear and the tidal field 
Ef^^, and the coupling between electric and magnetic 
Weyl curvature via the curl terms, which supports the 
propagation of gravitational waves. This set of equations 
would have been signficantly more complicated had we 
not chosen u^ to be comoving, with extra terms involving 
the momentum density and acceleration a^. In addition, 
the initial conditions must satisfy a set of initial-value 
constraint equations, although we will not need the ex- 
plicit form of the constraints here. 



B. 1 + 1 + 2 decomposition 

Under spherical symmetry, each comoving-orthogonal 
slice contains a preferred spacelike congruence with radial 
tangent vector r^, where r^r^ = 1. By analogy with the 
tensor ft.''^, defined in Eq. (1), which projects into the 
slices, we can define a tensor s^j^ by 



5/^ix — i^^v 



(16) 



which projects orthogonally to both u^ and r^. A 2- 
surface lying in a slice and orthogonal to r^ will be called 
a sheet. (The sheets are ordinary 2-spheres under spher- 
ical symmetry.) 

Any symmetric, trace-free tensor T^j/ which is orthog- 
onal to u'^ (i.e. which is a 3-tensor) can be decomposed 
according to 
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Here T = r^r^T^i, is a fully radially projected 2-scalar, 
T^ = sf^^rxT"^ is a projected 2-vector, and T^'' = 



^s^ry\-s^^'^s^x/2 



T 



kX 



^ „ OKA/^i -^ is a symmetric, trace-free 2- 
tensor. Any 3-vector can be similarly decomposed, al- 
though we will not need the explicit result. Any 3-scalar 
such as the energy density p is automatically also a 2- 
scalar. 



The utility of this 1 + 1 + 2 decomposition becomes 
immediately apparent when we realize that, under spher- 
ical symmetry, all 2-vectors and 2-tensors must vanish. 
Then, with straightforward manipulations, the evolution 
Eqs. (11) to (14) become 
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Here I have defined the 2-scalars S = 
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(18) 



(19) 



(20) 



(21) 



and 



E = r^^r'^Efj^i,. The set of tensorial partial differen- 
tial equations, (11) to (15), has been replaced with a 
simple, closed set of scalar ordinary differential equa- 
tions (ODEs). The magnetic Weyl tensor vanishes under 
spherical symmetry: no gravitational waves are possible 
in this case. 

An equivalent way to see that all quantities reduce to 
2-scalars in the spherically symmetric case is that any 
3-tensor or 3-vcctor can only be constructed from r^ and 
s^"; no other objects are available. [Eq. (16) relates h^i, 
to r^ and s^^.] Thus, e.g., the trace-free shear tensor cr^^ 
must take the form 



1 



'fJ.L' 



(22) 



C. Exact solution 

The set of ODEs, Eqs. (18) to (21), could be numeri- 
cally integrated, given initial conditions that satisfy the 
constraint equations. However, as is well known, an exact 
solution exists for the spherically symmetric dust space- 
time [8, 9, 10]. It is normally expressed in terms of the 



metric, but can be readily rewritten in terms of the co- 
variant 2-scalars p, 9, E, and 8. To express the exact 
solution, we must first introduce a time coordinate t that 
labels the comoving-orthogonal slices. Because the co- 
moving worldlines are geodesic, we can choose t to mea- 
sure proper time along the worldlines; henceforth this 
choice will always be assumed. Similarly, it is useful to 
define a radial coordinate, r, which is constant along the 
comoving worldlines. 
The result is 
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(23) 



(24) 



(25) 



The decaying mode corresponds to fluctuations in the 
bang time function, tsir), and can be set to zero by 
choosing tBir) = const. The growing mode corresponds 
to fluctuations in the spatial Ricci curvature '^'^^i? of the 
comoving slices, and can be made to vanish with the 
choice ^^'R ~ 0. It will be very important to eliminate 
the decaying mode in setting up realistic void profiles, 
since otherwise at early times the Universe becomes ex- 
tremely inhomogeneous, contradicting the standard pic- 
ture of extreme post-inflationary homogeneity. 

Although it will not be needed in this work, to facilitate 
comparison with previous studies the LTB metric can be 
written explicitly as 



ds^ 
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a 
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dr^ +Y^dn'^. 



(28) 



It reduces to the FRW metric in the homogeneous case 
(where the shear and electric Wcyl scalars, E and £, both 
vanish), and to the Schwarzschild metric outside some 
radius r = L, if p{r) = for r > L. 
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where Y — Y{r,t) is given implicitly by 



M 
Y = —{l-coshij), t-tB = 



M 



- (sinh r; — rf) , 



(27) 
and M = M{r), K = K{r) < 0, and ts = is(r) are ar- 
bitrary radial functions. Here the prime symbol denotes 
d/dr. (Related solutions are known for the cases K >0, 
but will not be needed here.) These expressions for the 
2-scalars were first written down in [24] . Given values of 
t and r for which a solution is desired, and choices for 
the three arbitrary functions M{r), K{r), and tB{r), we 
can solve Eq. (27) for Y{r,t) and its derivatives. Then 
we can evaluate each of Eqs. (23) to (26) to complete the 
solution. It is straightforward to verify explicitly that 
Eqs. (23) to (27) solve the DDEs (18) to (21) (as well as 
the initial- value constraints) for arbitrary Af(r), K{r), 
and tB{r). 

Because we are free to reparametrize the radial coordi- 
nate r, only two of the radial functions M{r), K{r), and 
tsir) describe physically distinct solutions. Therefore, 
the general spherically symmetric dust solution depends 
on two free radial functions. It will be crucial to un- 
derstand the physical significance of these two functions, 
which is the subject of the appendix, but we can see 
readily from the above exact solution that as f — > t^ (*") , 
for fixed r, the density p -^ oo, i.e. we approach a physi- 
cal singularity. This singularity is analogous to the initial 
singularity in the homogeneous and isotropic FRW space- 
time, but with extra radial dependence; hence the name 
"bang time" for the function tsir). 

In the appendix, it is shown that the two free radial 
LTB functions correspond to the growing and decaying 
mode for linear perturbations about FRW backgrounds. 



III. SCALAR PERTURBATIONS ABOUT LTB 
BACKGROUNDS 

A. Scalar-tensor coupling 

The problem of evolving linear perturbations on a 
spherically symmetric dust (LTB) background is in prin- 
ciple straightforward: we must linearize the set of evo- 
lution equations (11) to (15) about an arbitrary LTB 
background solution. In doing so, the symmetry of the 
background implies that any 2-vector or 2-tensor quanti- 
ties must be of first order, and hence their products can 
be ignored. 

First consider a 2-scalar subset of the full set. namely 
Eqs. (11) and (12) and the contraction between r^r" and 
Eqs. (13) and (14). Since we must relax (at linear order) 
the condition of spherical symmetry, the vector field r*^ 
can no longer be chosen to be exactly radial. Neverthe- 
less, we can suppose that r'^ departs from radial only at 
first order. We will see in Sec. IIIB that 2-scalar quanti- 
ties arc invariant under first-order variations in the field 
r^, so we have no need to fix the field at first order. 
Straightforward manipulations then lead to the 2-scalar 
linearized set: 
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Here W^ is the 2-vector part of the 1 + 1 + 2 decomposi- 
tion of i?^,y; recah Eq. (17). The symbol V represents 
the sheet-projected covariant derivative, defined in anal- 
ogy with Eq. (2), but with sheet projection tensors s^i^ 
replacing /i^^; i.e., it represents the covariant derivative 
in the "angular" directions. Finally, e^^ = e^^xr^. This 
set of equations agrees with the corresponding linearized 
subset of the full set derived in [18]. Again, a set of 
initial-value constraint equations is not shown. 

Remarkably, this set of 2-scalar equations is identical 
to the corresponding exact set under spherical symme- 
try, Eqs. (18) to (21), with the exception of the term 
e^i^V^T^'' which couples £ to the magnetic Weyl tensor. 
Physically, this term couples scalar perturbation modes 
to tensor modes, represented by TJ^i, [41]. While it is well 
known that scalars and tensors decouple when linearizing 
about FRW, such a coupling is expected to be a feature of 
the evolution on more general backgrounds, since a non- 
trivial background allows one to construct linear scalar 
quantities via, e.g., contractions between the background 
and tensor modes. (Recall as well that scalars couple 
to tensors at second order about FRW backgrounds, es- 
sentially since the second-order variables propagate in a 
nonhomogeneous first-order background.) 

Unfortunately, because of the scalar-tensor coupling, 
the hnearized set, (29) to (32), is not closed, and to com- 
plete the description of the dynamics we must include the 
2-vector part of the H^^ evolution equation, Eq. (15). 
This equation in turn couples to 2-vector and 2-tensor 
parts of the shear and electric and magnetic Weyl ten- 
sors, resulting in a very large set of equations [18] . On top 
of the added complexity of several more equations, the 
coupling to H^j^i, means that we must evolve a set of par- 
tial differential equations. If the coupling were absent, we 
would only need to solve a much simpler set of ordinary 
differential equations. Indeed, as I mentioned above, in 
that case the set of equations would be identical to the 
corresponding set under spherical symmetry, and hence 
we could employ the known exact solution to those equa- 
tions! (Of course, since we have ignored second-order 
terms in Eqs. (29) to (32), the use of the exact solution 
in this way could only be trusted up to first order.) 

The vast simplification of the dynamics in the absence 
of the tensor coupling leads naturally to the question: 
Under what conditions could the tensor coupling be ig- 
nored to a good approximation? It is certainly reason- 
able to consider the case where tensors are negligible at 
early times, before the LTB background becomes non- 
linear (e.g. at the time of last scattering), since many 
viable inflation models predict a very small primordial 
tensor contribution. However, even in this case, tensor 
modes will generally be sourced by scalars at late times, 
when the LTB background becomes nonlinear, according 
to Eq. (15). 

Nevertheless, in the case of interest, namely that of 
linear scalar perturbations on top of a mildly nonlinear 
spherical inhomogeneity, it is not expected that signifi- 
cant tensors will be sourced. It is expected that tensors 



would be sourced at times late enough, or scales small 
enough, that the scalar perturbations themselves have 
become nonlinear [25]. Essentially, nonlinear collapsing 
overdensities during structure formation are generically 
nonspherical, and hence source tensors through a signif- 
icant time-dependent quadrupole moment. Instead, in 
the case of interest in the present work, wc have linear 
fluctuations on top of a nonlinear background which itself 
cannot source tensor modes. 

In Ref. [26], the authors performed a numerical study 
of the linear stability of H^^ = in a few exact 
backgrounds, including a planar inhomogencous Szekeres 
model. They found that vanishing magnetic Weyl cur- 
vature was stable during pancake collapse, but unstable 
during collapse towards a spindle geometry. This sup- 
ports the claim that tensors can be ignored in the present 
context, where in the late stages of evolution of a void, 
the outer regions tend to collapse to form an ovcrdcnsc 
shell, which locally approximates pancake collapse geom- 
etry 

Note that dynamics with vanishing pressure and mag- 
netic Weyl curvature are generally described by a set of 
ODEs along the comoving worldlines, as can be seen di- 
rectly from the exact set, Eqs. (11) to (14). This was 
noted some time ago [27], and such a spacetime was 
dubbed silent. This name refers to the fact that without 
pressure and H^j^^, sound waves and gravitational waves 
cannot be supported, respectively, and hence no direct 
communication can exist between neighbouring world- 
lines. Apart from a few special cases, such spacetimes 
are thought to be generally inconsistent at the nonlin- 
ear level, in that the condition Hp_u = imposes a new 
constraint on the dynamics which is not consistent with 
the remaining dynamical equations [28]. However, in the 
present context wc are only assuming that the coupling 
to tensors can be ignored at linear order. Also, we do not 
require that the full tensor H^,^ vanish, only the weaker 
condition that e^^V^W can be ignored. 

Henceforth I will assume that the tensor coupling term 
can be ignored in Eq. (32) in evolving perturbations on 
LTB backgrounds. As a consistency check, I show in 
Section IV that for a particular LTB profile chosen to 
provide a rough fit to the luminosity distance-redshift re- 
lation of standard A cold dark matter (ACDM) models, 
the background shear, through which tensors must cou- 
ple to density perturbations, only very weakly effects the 
density perturbations. Such a test should be performed 
for any such calculation of scalar perturbations on LTB 
backgrounds. Obviously the full set of equations, includ- 
ing tensors, could in principle be evolved to check the 
validity of the weak tensor-scalar coupling approxima- 
tion. This might be necessary if we were interested in 
evolving perturbations on spherical overdensities at the 
late stages of collapse, instead of underdensities. 



B. Defining perturbations 

The dynamical equations derived in Section III A de- 
scribe the evolution of the total quantities p, 6, E, and 
£, i.e. backgrounds plus perturbations, in the linear ap- 
proximation. For many purposes it is very useful to ob- 
tain the evolution of appropriately defined perturbations 
alone. To do this in the metric-based approach to pertur- 
bation theory, one first provides a mapping (i.e. makes a 
gauge choice) between the real spacetime and a fictitious 
background (usually FRW). Perturbations in any quan- 
tity arc then defined as the difference between the exact 
value at some event and the background value mapped 
to that event. The freedom to vary the mapping results 
in the familiar gauge ambiguity. 

On the other hand, within the 1-1-3 covariant approach 
to perturbation theory (see [19] for a recent review) per- 
turbations are usually represented by spatial gradients 
orthogonal to the timelike direction u^^; e.g., density fluc- 
tuations are characterized by the 3-vector Df^p. The in- 
tention is to describe fluctuations in a gauge-invariant 
and covariant manner. However, it is easy to see that 
precisely the same ambiguity exists in this approach as 
in the metric approach. If wc choose w'' to be orthogo- 
nal to hypcrsurfaces of constant energy density, then we 
have Dfj^p = trivially, i.e. we have "gauged away" the 
density perturbation. Of course, as we have seen here, 
in the case of a single-component matter source, there 
is a natural choice for u^ which simplifies the dynamical 
equations, namely the comoving choice. But this choice 
is not necessary. 

While the ability to define perturbations entirely 
within the true spacetime is an advantage of the covari- 
ant approach for FRW backgrounds, it does not appear 
to be possible to extend this approach to LTB back- 
grounds. This is because, on the natural comoving- 
orthogonal slices, gradients such as I?^p do not generally 
vanish at background level, and hence they cannot be 
said to characterize the perturbations alone. Therefore, 
here I will employ a more standard approach, defining 
the density perturbation as the difference 



(5p(a;^) = p(x^)-/"'(r,<). 



(33) 



Here p{x^) is the density at event x^ in the true, per- 
turbed spacetime, and p^^-* (r, t) is the density at the cor- 
responding event in the background, in this case a spher- 
ically symmetric LTB spacetime. I define perturbations 
in the 2-scalars 0, S, and 8 analogously. The comoving 
worldlines in a dust spacetime provide a natural choice 
for the time coordinate. This coordinate is chosen (in 
both true and background spacetimcs) such that hyper- 
surfaces of constant t are orthogonal to the comoving 
worldlines, and t measures proper time along those world- 
lines. This fixes the "temporal gauge." 

There is inevitable freedom in choosing the radial coor- 
dinate r, since no radial vector r^ is naturally defined in 
the perturbed spacetime, and there is no unique notion 
of distance in a curved spacetime. We can defeat this 



ambiguity by considering the intended application of the 
perturbed LTB formalism, namely the evolution of linear 
scalar perturbations on top of a spherical void. Such a 
void becomes nonlinear at late times, but, as discussed 
in the appendix, must consist of a pure growing mode, to 
be consistent with the inflationary paradigm. Hence the 
spacetime approaches FRW plus linear perturbations at 
early times (such as the time of last scattering). Wc can 
define r in the FRW background at such early times to 
give radial proper distance from some centre. Then we 
can define r on the perturbed spacetime at early times 
to be any radial coordinate that agrees, to first order, 
with the background radial coordinate. Taking the back- 
ground at this early time to be homogeneous FRW, any 
linear variations r ^> r -\- 5r in the radial coordinate can- 
not affect the values of any perturbed physical quantities, 
to first order. [Physical perturbations in FRW are invari- 
ant under purely spatial gauge transformations, at first 
order [29, 30].) We then define r at later events such that 
it is constant along each comoving worldline in the true 
spacetime, i.e. r is "dragged" along by the fluid flow. 
(Therefore comoving worldlines in the background are 
mapped to comoving worldlines in the perturbed space- 
time.) The spacelike vector field r^ is defined to be nor- 
mal to surfaces of constant r. These surfaces will be 
linearly perturbed versions of 2-spheres. 

Note that a type of radial gauge dependence may be 
present in perturbations at late times. In regions where 
the background energy dcnstity, e.g., has a significant 
spatial gradient, it may be possible to gauge away a den- 
sity perturbation 6p with a radial gauge transformation. 
(Indeed this gauge dependence may continue to early 
times if we consider the background to be LTB, rather 
than FRW, at the early times.) However, we will not be 
able to gauge away 5p everywhere, since the density gra- 
dient will not be large everywhere {dp/ dr must vanish 
at the origin if p is to be smooth there). In the FRW 
case, the usual monotonicity of the background p(t) does 
allow us to set 5p ^ Q everywhere. The possible gauge 
dependence in the LTB case must be kept in mind when 
interpreting the results of calculations. Of course, a cal- 
culation of any observable quantity must be independent 
of any gauge choice. [An example would be the variance 
of 5p across a 2-sphere of constant redshift (which is di- 
rectly observable), rather than some arbitrary coordinate 
r (which is not).] 

Recall that in the 1 + 1-1-2 formalism, 2-scalars are 
constructed from 3- vectors or 3-tensors by fully project- 
ing along r^. Thus we might be concerned that linear 
variations in r^, which will necessarily follow from lin- 
ear variations in the coordinate r, might change the val- 
ues of 2-scalars, at least at late times when the back- 
ground has become nonlinear. However, this does not 
happen, at first order, because of the isotropy of the back- 
ground. Consider an exact, spherically symmetric LTB 
background, with exactly radial vector r^ and linearly 
perturbed vector 



(5r^. 



(34) 



Imposing the normalization ff^ = 1 gives f r^ = 1 + 
0(2). Therefore, defining the shear scalar E using the 



radial ?'^. S 



r^r^cFfi^y, we have 



ft^r^r 



0{2). 



(35) 



(A similar result applies to the electric Weyl tensor.) 
That is, linear variations in r^ do not affect the values of 
2-scalars, to first order. Clearly the presence of pertur- 
bations will not change this result, since the variations 
in the values of 2-scalars due to the presence of linear 
perturbations as r^ is varied at first order will also be of 
at least second order. 

To summarize, the temporal gauge is chosen to be co- 
moving, while the radial gauge is any linear variation 
about "initial slice proper distance gauge." The tempo- 
ral choice simplifies the dynamical equations, while the 
radial choice enables the simple adaptation of standard 
results from the theory of perturbations on FRW back- 
grounds, as we will see in Sec. HIE. 



C. LTB transfer functions 

Inserting the definition of 5p, Eq. (33), and the cor- 
responding expressions for the other 2-scalar perturba- 
tions, into the linearized set, Eqs. (29) to (32) (dropping 
the tensor coupling), and discarding terms quadratic in 
the perturbations, we obtain 



8p = -eSp 



5Y. 



AirGSp - 3E(5E, 



E ST. 



5£ 



T.\5£-£ 



-T6e - 58, 



-5Y. 



(36) 



(37) 



(38) 



47rG'(/9(5E + E(5p). 



(39) 



In obtaining these equations, the fact that comoving 
worldlines in the background are mapped to comov- 
ing worldlines in the perturbed spacetime has allowed 
us to subtract off the background equations of motion, 
Eqs. (18) to (21). Also, note that it is valid to trivially 
perturb time derivatives (e.g., 5p = p — p^^'), since t 
measures proper time in both background and perturbed 
spacetimes. 

In a sense, this completes the problem of perturbations 
on LTB backgrounds, since the set of ODEs, Eqs. (36) 
to (39) , can be readily solved numerically (or even using 
the exact solution from Sec. II C, remembering that that 
solution can only be trusted to linear order here). How- 
ever, with applications in mind, it will be very useful to 



develop a description of the evolution of the perturba- 
tions in terms of a new set of transfer functions. This 
will allow us to derive results that parallel remarkably 
closely the standard results from perturbation theory on 
FRW. 

Consider the column vector of 2- vector perturbations, 
dXi{t,r,n^), where Xi = p, 9, E, £. Here n^ is the 
spatial unit vector from the origin, r ~ 0, towards the 
point at which 6Xi is evaluated, so (r,n'^) are the polar 
coordinates of that point. Then the linear perturbations 
at some arbitrary time t must be related to those at some 
initial time ti via the matrix product 



5X,{t, r, n'^) = T„(f, r)5Xj{t„r, n^). 



(40) 



where Ty (t, r) is a 4 x 4 matrix of transfer functions. Lin- 
earity of evolution demands that Eq. (40) take a linear 
form. The property of "silence," i.e. that the perturba- 
tions are solutions of ODEs, means that the solution at 
event (t, r, n^) can only depend on the initial condition 
on the same worldline at (ii, r, n^). (Initial velocities 5Xi 
are not needed since the ODEs are first order.) Finally, 
the isotropy of the background implies that the transfer 
functions T^ must depend only on t and r. 

In the case of perturbations on FRW, the correspond- 
ing transfer functions are generally scale (i.e. k) depen- 
dent, although this is only true when pressure is present 
to support acoustic waves and introduce spatial deriva- 
tives into the dynamical equations. In the LTB case here, 
no such scale dependence exists. Also, the FRW transfer 
functions appear to be considerably simpler than their 
LTB counterparts, in that they depend on t alone in the 
homogeneous case. However, again with applications in 
mind, we will very often be interested in evaluating the 
perturbations on the observer's past light cone (with the 
observer located at the centre of symmetry). The sur- 
face of that cone defines a definite relationship r = r{t) 
(at background order), and so the LTB transfer functions 
can in practice be considered to be functions of the time 
alone, r,;j(i,r) =T,j{t,r{t)). 

Thus the problem of linear perturbations on LTB back- 
grounds, which at first may have appeared quite daunt- 
ing, is reduced to that of calculating a small number of 
transfer functions of a single variable, for most practi- 
cal applications. As we will see next, the problem can 
be simplified even further by consideration of the initial 
conditions. 



D. Initial conditions 

In order to use Eq. (40) to evolve the perturbations 
we must first specify initial conditions. To do this, first 
recall again the context in which we are considering LTB 
backgrounds here. We require growing mode solutions, 
which approach linearly perturbed FRW at early times. 
Therefore, we can set up initial conditions for the fluc- 
tuations at some sufficiently early time ti, exactly as in 



the standard FRW case. For definiteness (and well sup- 
ported by current observations [31]), I will assume that 
the primordial fluctuations are scalar and adiabatic. In 
addition, with their inflationary origin in mind, I will take 
the perturbations at the initial time, U, to be purely in 
the growing mode (this assumption is also supported by 
the observations [32]). Since the above calculations have 
assumed vanishing pressure, we must take the initial time 
for the specification of perturbations to be sufficiently 
late that matter domination is a good approximation. 

With these assumptions, we can reduce the number of 
LTB transfer functions that need to be calculated, since 
some of the initial perturbations can be written in terms 
of the others. It will be particularly useful to write all 
of the initial perturbations in terms of the comoving cur- 
vature perturbation, TZ{x^). This quantity is constant 
on super-Hubble scales (for adiabatic modes), and hence 
is very convenient for comparing with the predictions of 
a particular inflationary model. Constraints on the pri- 
mordial value of TZ are also explicitly provided in obser- 
vational results (see, e.g., [31]). 

During matter domination, the comoving energy den- 
sity perturbation can be written in terms of TZ as 



5p_ 
P 



18 D' 



n 



(41) 



(see, e.g., [33]), where D'^ = D^D^ is the spatial Lapla- 
cian. At the initial time <i, the shear tensor for scalar 
fluctuations can be derived from a scalar function cr(a;''), 
according to 



CT^i/ = D^D^a - -h^^D'^a. 



(42) 



Again during matter domination, the comoving shear 
scalar a can be written in terms of the comoving cur- 
vature perturbation as 



Qn 



(43) 



(see, e.g., [30]). Combining these results we can write for 
the shear 2-scalar initial condition at i^, 

e ho-^ydr-' 3 y ^ ' 

Note that, even though both 5p and 5Y. are completely 
determined by the comoving curvature TZ, they are not 
simply proportional. Geometrically, it is possible, e.g., 
that some localized linear fluctuation in TZ produces a 
5p at some point, but no (5E, if the fluctuation is sym- 
metric about the point. Nevertheless, in the appendix, 
Eqs. (A.l) and (A. 2) show that the growing modes for 
the comoving gauge energy density and expansion per- 
turbations are related by 



P 



(45) 



Similarly, using Eq. (38) (for the relevant case at ti that 
the background shear S is first order and hence can be 
dropped from this equation) and Eq. (44) , we can relate 
the electric Weyl 2-scalar to the shear via 



5£ 



1 



-MH. 



(46) 



It is important to point out that these initial conditions 
must automatically satisfy the Einstein initial-value con- 
straint equations, to linear order. Hence all that remains 
is to determine the transfer functions to evolve the initial 
conditions. 

The proportionalities (45) and (46) apply at the initial 
time, ii, and hence they reduce the dimensionality of 
the matrix of transfer functions in Eq. (40) to 2 x 2. 
For example, the evolution of the density perturbation is 
given by 






Pit,r) 



Tp^it.r) 



P{U) 



(47) 



and we only need to calculate two transfer functions. In 
this expression I have scaled the perturbations so that 
the transfer functions are dimensionless and positive, and 
the function Tpp approaches unity at early times {Tp^ 
approaches zero at early times). 



E. Harmonic decomposition 

One of the difficulties that comes to mind when first 
considering the problem of linear perturbations on LTB 
backgrounds is that of harmonic expansion. A crucially 
important tool in perturbation theory on flat FRW back- 
grounds is the ability to expand fluctuations in Fourier 
modes or spherical harmonics and spherical Besscl func- 
tions. This approach can be generalized to the case of 
spatially curved FRW backgrounds. But an LTB back- 
ground can have an arbitrary radial dependence of the 
spatial curvature, so it appears hopeless to look for ra- 
dial harmonics generally (although the standard spheri- 
cal harmonic functions, the Yfm's, can be used to expand 
the angular dependence, of course). 

Again, though, we are considering spacetimes which 
approach linearly perturbed FRW at early times. There- 
fore we should be able to employ the usual harmonic ex- 
pansions at a sufficiently early time. Then we can evolve 
the perturbations using the transfer functions introduced 
in Section IIIC. The use of standard harmonic expan- 
sions at early times will, e.g., allow us to answer statisti- 
cal questions about the perturbations at late times, using 
the standard formalism of Gaussian random fields, just 
as we do on FRW backgrounds. 

To see how this works explicitly, consider first the co- 
moving curvature perturbation TZ. At the initial time ti , 
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when the fluctuations are well described by linear pertur- 
bations from FRW, we can expand TZ in spherical har- 
monics as 



where jt is a spherical Bessel function of the first kind 
and k is a, comoving wave number. The curvature pertur- 
bation TZ evaluated at time ti, during matter domination, 
is related to the primordial, constant value HP'' via the 
standard matter transfer function T{k): 

7^ft„(fc,t,) = T(fc)7^P'„(fc) (49) 

(see, e.g., [33]). The function T[k) describes the suppres- 



sion of power in TZ incurred while each mode is inside the 
Hubble radius during radiation domination. Assuming 
Gaussian random primordial fluctuations, the statistics 
of T^P'' are entirely encoded in the standard expression 



(48) {nZ,{k)TlT^,{k')) 



2^.Mfc) 

fc3 



5{k-k')5u,5,nm'. (50) 



Here the angled brackets indicate averaging over the en- 
semble of realizations of the fluctuations. For a scale 
invariant primordial spectrum the dimensionless power 
spectrum T'-iz{k) is constant. 

Taking the energy density as an example, by combining 
Eqs. (41), (44), and (47) to (49), we can write 



5p{t,r,n^') 2 [2 



p{t.r) 



J- fdk kTik)Y,'R'Tmik)Tp{t,r,k,e)Y,minn, 



where I have defined a new transfer function by 



Tp{t,r,kJ) 



QiHi 



Tpp{t,r)Mkr) - Tpj:{t,r) ( j'^ikr) + -je{kr) 



(51) 



(52) 



Here I have explicitly included the scale factor ai at 
time ti, and Hi = 0(ti)/3. The notation j" indicates 
the second derivative oi ji with respect to its argument. 
Note that the details of a particular void profile enter 
these expressions only through the LTB transfer func- 
tions Tpp{t,r) and Tps{t,r). I stress that the decompo- 
sition into wave numbers k in these expressions has been 
done at the early time ti, so the harmonic expansion is 
valid. 

We can now calculate any statistical aspect of the 5p 
field, at any time. For example, consider the angular 
power spectrum of density fluctuations on a sphere of 
some particular radius r and at some particular time t, 
defined by 



Cf''{t,r)^ 



5p{t,r,n^) 
P{t,r) 



Y.l.inndnj 



(53) 



[presumably r ~ r{t) if we are interested in observations 
on the past light cone]. Using Eqs. (50) and (51), and 
the orthonormality property of the spherical harmonics, 
we can evaluate this expression as 



Cf'it^r) 



167r 
'25 



dk 



T\k)T^Jt,r,kJ)Vn{k). (54) 



Using similar techniques we can calculate a large variety 
of observable quantities, which allows us to confront void 
models of acceleration with observations at the level of 
perturbations [21]. 



IV. NUMERICAL EVOLUTION 

It is now a straightforward exercise to calculate any 
of the LTB transfer functions defined above. First, we 
must specify an LTB background about which to perturb. 
For numerical convenience I have used the function K{r) 
defined in Sec. II C to specify the intial profile, using 
the exact solution, Eqs. (23) to (27). As discussed in 
the appendix, I have set the bang time function tsir) 
to zero to eliminate the decaying mode. The remaining 
radial gauge freedom was used to set M{r) = r^, which is 
the value the corresponding function would take in FRW; 
although as discussed in Sec. Ill B any choice of radial 
coordinate that agrees with the background at zeroth 
order would do. 

For the initial profile I chose 



K{r) 




r<L, 



r> L. 



(55) 

This profile is characterized by two parameters, a width 
L and an amplitude Km- For r 3> i, it approaches flat, 
matter-dominated FRW, i.e. the Einstein-de Sitter (EdS) 
universe. I chose units such that the initial time was 
ti ~ I, and an amplitude which gave a spatial curvature 
of '^^^R/9^ — —5.5 X 10^"* at t,, so that the initial profile is 
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a small perturbation from FRW. The profile was chosen 
to roughly fit the luminosity distance-redshift relation 
for standard ACDM models, while providing a reason- 
able density parameter £!,„ today. However, evaluating 
the goodness of fit of the current LTB model to data is 
beyond the scope of this paper, and will be the subject 
of a future detailed study [21]. 

In Fig. 1, I display the background profiles for the en- 
ergy density, expansion rate, 2-scalar shear, and spatial 
curvature '■^■'i? calculated from the linearized Einstein en- 
ergy constraint equation. 



^l^SnGp-l^'^R- 



(56) 



All quantities are displayed at the late time t = 1.0 x 10^, 
and were calculated using the exact solution, Eqs. (23) 
to (27). In the plot, the expansion and shear 2-scalar 
are normalized by the value of the expansion rate at the 
origin, 0,nax- The energy density and spatial curvature 
are normalized such that their difference should be unity, 
according to Eq. (56) (the shear squared term is small 
compared to the density and curvature). 
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Figure 2 shows the transfer functions Tpp and T^s, 
which determine the evolution of the density perturba- 
tion according to Eq. (47). They were calculated by nu- 
merically evolving the linearized set, Eqs. (36) to (39). 
They have been normalized via 



Tinorm 
PP 



it,r) 



2/3 



j] Tpp{t,r), 



(57) 



and similarly for Tp^. This normalization was chosen so 
that, for the EdS case, T^°™it) ^ 1 [and TJ^i"^{t) = 0] 
for all time. The transfer functions are shown at the 
logarithmically-spaced times t = 46.4, 2154, and 1.0 x 
10^ 




FIG. 2: Transfer functions Tpp and T^e, normalized so that 
ynorm _ -^ ^^ g^jj tlmes foi an Einstein-de Sitter background. 
The functions are presented for the times t = 46.4 (red, dotted 
line), t = 2154 (blue, dashed line), and t = 1.0 x 10^ (black, 
solid line). Significant suppression of perturbations occurs at 
late times within the void, and density perturbations are also 
sourced by shear fluctuations at late times through Tps. 



FIG. 1: Background LTB radial proflle at the late time t = 
1.0 X 10^. Energy density p, volume expansion rate 9, shear 2- 
scalar E, and spatial curvature ^^'R are scaled as described in 
the text. This profile corresponds to an underdensity near the 
origin with effective density and curvature parameters Qm ~ 
0.3 and Qk = 0.7. 



Note that the normalized energy density and spatial 
curvature plotted in Fig. 1 correspond directly to the 
FRW definitions of the matter and (negative of the) cur- 
vature parameters, Qm and fix- Therefore we can see 
that this LTB profile corresponds, at the late time plot- 
ted, to a significant underdensity extending from the cen- 
tre to r ~ L, with effective density and curvature param- 
eters fim = 0.3 and Hk = 0.7 at the centre. Correspond- 
ing to this underdensity is a higher expansion rate than 
the asyinptotic flat region. The shear is small except at 
the boundary region r '^ L; it must vanish at the centre 
and in the asymptotic FRW region by symmetry. 



We can see from Fig. 2 that at early times, and for 
r ^ L, the transfer functions approach their EdS values. 
That is, the density perturbation is evolving as expected 
when the background is (locally) nearly EdS. However, 
we see a significant suppression in the growth of density 
perturbations at late times within the void. This sup- 
pression has an intuitive origin: near the centre of the 
void, the spacetime approximates the open FRW geom- 
etry, for which it is known that the growth of structure 
is suppressed compared with EdS (see, e.g., [33]). It is 
possible to verify that the amount of suppression at the 
origin quantitatively matches the open FRW prediction 
[21]. 

The function Tp^: plotted in Fig. 2 illustrates that, 
at late times in the nonlinear background regime, shear 
fluctuations can source density fluctuations, via coupling 
through the background shear [recall Eqs. (36) and (37)]. 
The importance of this effect is not clear from this fig- 
ure, however, since the two transfer functions Tpp and 
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Tp-£^ multiply different combinations of Besscl functions 
in Eq. (52). 

Finally, Fig. 3 illustrates the transfer function 
Tp{t,r,k,£), which was defined in Eq. (52), again nor- 
malized according to 

\t,r,k,£) 



rpnormf 
P 



2/3 

^) Tpit,r,k,t), (58) 



for the case ^ = 10 (the general features are very sim- 
ilar for other values of £). With this normalization, we 
have Tp{t, r, k, £) = [k / aiHi)'^ ji{kr) for all t and r for the 
special case of the EdS universe, and this case is plot- 
ted in Fig. 3. Also plotted arc the curves for the LTB 
model defined above, and evaluated within the periphery 
of the void, at r/L ~ 0.7, and near the void centre, at 
r/L = 0.01 (cf. Fig. 2). For the LTB cases, the transfer 
function is evaluated on the past light cone. Figure 3 
illustrates again the suppression of matter fluctuations 
near the origin with respect to the EdS case. 
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FIG. 3: Normalized transfer function T^°'-'"'{t,r,k,£) for £ = 
10. The curves are shown for a point within the periphery of 
the LTB void, at r/L = 0.7, as well as near the void centre 
at r/L = 0.01. Also shown is the local FRW approximation 
(LFRW) for r/L = 0.7, together with the Einstein-de Sitter 
case (EdS). 



To illustrate the effect of the background shear E and 
tidal field £, which are non-negligible only within the 
peripheral region of the void (recall Fig. 1), Fig. 3 also 
presents the transfer function Tp(t,r,k,£) evaluated at 
r/L = 0.7, but calculated by ignoring the background 
shear and tidal field in the linearized evolution Eqs. (36) 
to (39). This approximation can be called the "local 
FRW approximation," since setting S = f = re- 
sults in the evolution of the local quantities (background 
and perturbations) along each worldline that is identical 
to that of some FRW universe. We see that the local 
FRW approximation is very good, even in the void pe- 
riphery, where the background shear and tidal field are 
largest. (The local FRW approximation transfer function 
for r/L = 0.01 is indistinguishable from the correspond- 
ing curve plotted in Fig. 3.) Note that in the local FRW 



approximation we have Tps = 0, so Fig. 3 shows that 
modifications to the matter power spectrum due to shear 
fluctuations (as encoded in the transfer function Tp^) are 
subdominant to the main suppression which is captured 
well by the local FRW approximation. 

Recall that the tensor modes, which we have ignored 
in all of these calculations, couple to Sp only through the 
shear [see Eqs. (29) to (32)]. Since the magnitude of the 
background shear does not depend on the amplitude of 
linear tensor perturbations, the weak effect of the shear 
that we have observed (i.e. the excellent accuracy of the 
local FRW approximation) means that it is very unlikely 
that the tensors would have a significant effect on the 
density fluctuations had the tensors been included. It 
should be stressed that this result will depend on the 
background profile: a background more strongly nonlin- 
ear, and hence with larger shear, than the one studied 
here will couple both scalar shear fluctuations and ten- 
sor modes nrorc strongly to the density perturbations. 
Thus the importance of the background shear should be 
evaluated whenever evolving scalar perturbations on an 
LTB background. As well, the effect of the background 
shear may be larger on quantities other than the density 
perturbations. 



V. DISCUSSION 

It is expected that the formalism developed here could 
be applied to constraining void models for acceleration 
with a wide variety of structure-related data. Such a 
program could perhaps be dubbed "testing homogeneity 
with inhomogeneities." Examples include the shape (and 
amplitude) of the local matter power spectrum, a rigor- 
ous treatment of baryon acoustic oscillations, and the ef- 
fects of gravitational lensing by linear structures. These 
applications are currently under study [21]. Apart from 
void models, it is expected that the present approach 
could flnd use elsewhere, in particular in examining the 
early stages of the standard structure formation process. 

A particularly interesting application would be a full 
calculation of the effects of a void inhomogencity on the 
CMB. The effect of varying the angular diameter dis- 
tance to the last scattering surface in standard FRW cos- 
mologies is a change in the angular peak scale (indeed, 
this is the dominant effect on the spectrum from vary- 
ing the time of observation, keeping everything else fixed 
[34]), which can already be calculated at background level 
in LTB models [11]. The free-streamed primary CMB 
anisotropy power spectrum is expected to be very robust 
to the presence of isotropic inhomogencity. However, re- 
call from Fig. 2 that the presence of a large nonlinear void 
results in a suppression of the power of scalar fluctuations 
inside the void. Such a suppression exists at late times 
in standard ACDM models (as well as in open models), 
where it leads to the integrated Sachs- Wolfe (ISW) effect. 
An analogous effect should occur in void models, which 
will provide another handle on constraining these models. 
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The ISW effect modifies the largest angular scales, and 
so its presence is difheult to distinguish directly due to 
cosmic variance. However, correlations with large-scale 
structure have been measured with high confidence [35], 
and this is a signature that the present formalism could 
be used to calculate for void models. Even if the overall 
amplitude of the ISW-structurc correlations turns out to 
be insufficiently constraining, it would be very surprizing 
if the redshift-dependence of the effect would be similar 
between void and ACDM models. 

It is important to stress again that the present ap- 
proach is grounded within the well-developed theory of 
linear perturbations on FRW backgrounds, and many of 
the familiar standard techniques from that field can be 
applied here. The main difference is the presence of a 
small number of new, simple to calculate linear trans- 
fer functions. In practical terms, these transfer functions 
only need to be calculated once for a particular void pro- 
file, and can be expressed as a function of a single vari- 
able if results are required on the past light cone, as is 
normally the case. This connection with standard FRW 
perturbation theory was made possible by the crucial as- 
sumption that voids contain no decaying modes, so that 
they are consistent with the early homogeneity predicted 
by models of inflation. 

On the theoretical side, future work could include the 
calculation of the remaining 2-vector and 2-tensor shear 
and electric Weyl degrees of freedom, which, in the ab- 
sence of coupling to iJ^jy, also evolve according to ODEs, 
as Eqs. (13) and (14) show. Apart from possible obser- 
vational consequences of these components, this would 
enable a consistency check on the accuracy of ignoring 
the tensor coupling by monitoring the irrotational con- 
straint equation [19] 



Hf^i, = curlo-^,^. 



(59) 



However, it is important to stress that, for the void 
model studied here, the effect of the background shear 
on the density fluctuations was found to be weak. Thus 
any tensors produced arc not expected to have a signifi- 
cant effect on the density perturbations. The dominant 
effect, namely a significant suppression of power within 
the void, was found to be captured well by the local FRW 
approximation, which ignores the background shear and 
tidal field. Crucially, this dominant suppression near the 
origin is independent of the details of the void profile 
away from the origin, since the background shear must 
vanish at the origin by isotropy, so that the local FRW 
approximation is exact there. Similarly, the degree of 
suppression near the centre will not be affected by ten- 
sors, which require background shear to couple to density 
fluctuations. 

To conclude, proponents for void scenarios argue that 
they are conservative in that they require no new, myste- 
rious dark energy component (although, of course, they 
leave the solution of the cosmological constant problem 
to future work!), and that they may naturally explain 
the "coincidence problem," since structure started to go 



nonlinear only recently. Ultimately, though, wc will need 
to resort to the widest range of observations possible to 
determine if the apparently natural philosophical stance 
of the Copernican principle actually reflects our place in 
the Universe. 
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APPENDIX: GROWING AND DECAYING 
MODES 

1. Linearized FRW solution 

We saw in Sec. II C that the general spherically sym- 
metric dust (LTB) spacetime is specified by two free ra- 
dial functions. To help elucidate the physical significance 
of these two free functions, consider an LTB spacetime 
which is very "close" to FRW, in the sense that it can be 
treated as a linear perturbation from FRW. It is well 
known from the theory of cosmological perturbations 
(see, e.g., [36]) that a general linear dust perturbation 
on FRW is described by two modes, one growing and 
one decaying. [This can be seen easily from the exact 
Eqs. (11) and (12): For an FRW background, the shear 
term a^^a^^ is second order, and hence the linearized 
Eqs. (11) and (12) decouple from the remaining evolu- 
tion equations and form a second order system with two 
independent solutions.] Therefore we can conclude that 
the two free functions of an LTB spacetime correspond to 
growing and decaying modes in the linear regime about 
FRW. This point was made some time ago in Ref. [37] 
(see also [38]), but it will be useful to derive the result 
here in a more explicit manner. 

The importance of this result is that the presence of 
a significant decaying mode at late times will result in 
extreme inhomogeneity at early times, contradicting the 
standard inflationary scenario, according to which the 
Universe is expected to be extremely homogeneous at the 
end of inflation. Hence we will need to know how to set 
the decaying mode to zero in setting up a spherical LTB 
profile. This reduces the freedom available in choosing a 
profile, and hence should make it more difficult to fit all 
of the data. 

The evolution of the comoving gauge dust energy den- 
sity perturbation 5p in a pressureless FRW background 
is given by 



7^(^)"^-'^) + 7^*'- 



(AJ) 



where Di{x) and D2{x) are the amplitudes of the growing 
and decaying modes, respectively, and ig is an arbitrary 
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reference time. [It is straightforward to show that this 
expression is the solution to the Hnearized Eqs. (11) and 
(12).] Similarly, it can be shown that the comoving gauge 
expansion perturbation S0 [42] is given by a growing and 
decaying mode, 







2/3 



D^ix) + l^-fD2ix). (A.2) 



Using these two expressions and the linearized energy 
constraint equation, we can show that the 3-Ricci curva- 
ture of the comoving slices, '-^^R, evolves according to 



2 6*2 



2/3 



Di{x). 



(A.3) 



Note that the curvature perturbation consists of just the 
growing mode. The normalization of the 3-curvaturc per- 
turbation in Eq. (A.3) has been chosen to agree with the 
conventional FRW curvature parameter, 



Qk^- 



3(^ 
2 612 



(A.4) 



in the homogeneous limit. 



2. Vanishing growing mode 

I showed above that a comoving 3-curvaturc perturba- 
tion corresponds to a pure growing mode in linear the- 
ory on an FRW background. Therefore we can isolate 
the decaying mode by considering solutions with vanish- 
ing spatial Ricci curvature ^^^R. For an LTB spacetime 
which is a small perturbation from FRW, the exact en- 
ergy conservation and Raychaudhuri equations, (11) and 
(12), become, to first order. 



P 



-—-inGp, 



(A.5) 



(A.6) 



since the shear cr^j/ is a first-order quantity. Under 
the same approximation, the Einstein energy constraint 
equation becomes 



92 1 



(A.7) 



For the case '^^i? = 0, it is simple to verify that the 
solution to these equations is 



p(r,i) = 



6^G[t-iB(r)] 



2 ' 



(A., 



Here the function tB{r) has an arbitrary radial profile, 
but in order that the shear (and electric Wcyl tensor) be 
small, the size of radial fluctuations in Ib (?') must not be 
too large. 

To quantify this, consider a homogeneous background 
solution 



P^'\t) 



1 



GttG [t-t 



.m 



9^°\t) = 



t-t 



(0)' 



(A.IO) 



(A.ll) 



near the perturbed solution (A. 8) and (A. 9), where tg 
is a constant. Define the energy density perturbation by 



Ap(r,i)=p(r,t)-p(")(t). 



(A.12) 



Then, expanding in powers of Ats/t, where At sir) — 
tsir) —tg, we have 

^^-2^^ + Or^V, (A.13) 



n(0) 



to t 



t 



with to a-n arbitrary reference time. Comparing this ex- 
pression with Eq. (A.l), we see that we have indeed iso- 
lated the decaying mode, and its amplitude can be read 
off of Eq. (A.13) as 



A similar calculation gives 



,AtsW 
to 



Ae _ Atstp .^[Ati 

em to t [ t 



(A.14) 



(A.15) 



which, with Eq. (A.2), again gives Eq. (A.14). Therefore 
it is fluctuations of the bang time tsir) that give the 
decaying mode, and we can set the decaying mode to 
zero by setting i_B(r) = const. 

It is very important to point out that the comparison 
between the ordinary FRW perturbations Sp and SO and 
the perturbations Ap and A9 is valid, since the FRW 
perturbations were specified in comoving gauge, and the 
LTB solutions Eq. (A. 8) and (A. 9) are specified with re- 
spect to the comoving congruence. 

Finally, note that at times early enough that t ~ 
At sir), the linear approximation breaks down. This 
is expected, since the decaying mode diverges at early 
times. This result is intuitively clear, as we expect dif- 
ferences in "bang time" to become all-important at early 
times, but to be inconsequential at sufficiently late times. 



3. Vanishing decaying mode 



0(r,t) 



t-teir)' 



(A.9) 



Now that we have isolated the LTB decaying mode, we 
can examine the growing mode, by looking at solutions 
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with tsir) = 0, but with nonvanishing curvature of the 
comoving sHces. For definiteness I will consider the open 
case, ^^■'i? < 0; the other cases follow similarly. 
In this case, Eqs. (A. 5) to (A. 7) have the solution 



p(r,i) = po(r) (^— j , 



(A.16) 



Combining this with the lowest-order expression 



gives 



1-"^.O„3^0(^5) 



m 



3/2 
K.Si 



(A.23) 



0(r,i)=3-, 
a 



(A.17) 



where the quantity a = a{r,t) is given parametrically by 



a 
ao 



1-n 



K,0 



2ilK,0 



(cosh?/ — 1), 






(A.lSa) 



(A.18b) 



Quantities with subscript are evaluated at reference 
time < = to, and are related by 



{l-nK.fl)Ol^2AnGpo, 



(A.19) 



and. by analogy with the homogeneous FRW case, I de- 
fine 






(A.20) 



Therefore the solution is determined by one free radial 
function, which describes the spatial Ricci curvature on 
some reference comoving slice. 

Again, consider a spatially flat homogeneous back- 
ground solution 



p(°)(0 



6^Gi2 ' 



(A.21) 



and define the density perturbation Ap by Eq. (A. 12). 
If we are to evaluate this perturbation at time to, then 
for the linear approximation to be valid we require 
^Kfiir) «C 1, i.e. the perturbation must correspond to 
a small curvature fluctuation. Then Eq. (A. 18b) tells us 
that r] <^ 1 (this is equivalent to the "small u" approxi- 
mation in [39]). We can now solve the set Eq. (A. 18) for 
a/ao perturbatively in 77. Keeping the next-to-leading or- 
der terms in the expansions for the hyperbolic functions, 
I find 



ao 



— [dirGpot 



2n1/3 



1 1 



80 nKn 



^^'°r,'i + 0{if). (A.22) 



ao 



= {enGpof 



,1/3 



^K,0 



5 (1 - f]x,o)2/3 \t 



2/3' 



Combining this expression with Eq. (A.16) gives 



(A.24) 



Pir,t) 



1 



67rGt2 



1-^- 



^K,n 



2/3' 



5 (1 - f)K,o)^/3 \to. 

(A.25) 
Therefore the density perturbation is finally determined 
to be 



Ap 

P' 



(0)- -^"^^■"U 



2/3 



(A.26) 



at lowest order. Comparing with Eq. (A.l), we see that, 
as expected, Eq. (A.26) consists of a pure growing mode. 
As with the decaying mode case, we can read off the 
growing mode amplitude as 



5 



(A.27) 



With Eq. (A.20), we recover Eq. (A. 3), i.e. the growing 
mode is given by the spatial curvature, as expected. 

A similar calculation gives 



A6' 



l""""" [To 



2/3 



(A.28) 



which, with Eq. (A.27), is again consistent with the FRW 
expression for the growing expansion perturbation mode, 
Eq. (A. 2). Therefore, the growing mode for LTB space- 
times near FRW is indeed given by the spatial curvature 
of the comoving slices. 
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